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Abstract 
o 

^^ ' Heavy colored scalar particles, which exist in many models of new physics, can be pair produced 

^1 ■ at the Large Hadron Collider (LHC) via gluon-gluon fusion and possibly form quarkonium-like 

bound states. If the scalars are also charged under the electroweak gauge group, these bound 

Qh' states can then decay into electroweak bosons. This yields a resonant cross section for final states 

^' such as 77 that can exceed Standard Model backgrounds. This paper studies this process in 
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the Manohar-Wise model of color-octet scalars (COS). Important threshold logarithms and final 
state Coulomb-like QCD interactions are resummed using effective field theory. We compute the 
resummed cross section for gluon-gluon fusion to COS pairs at the LHC as well as the resonant 
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f"^ . cross section for octetonium decaying to 77. The latter cross section exceeds the Standard Model 



di-photon cross section when the COS mass is less than 500 (350) GeV for y/s = 14 (7) TeV. 
Nonobservation of resonances below these energies can significantly improve existing bounds on 
COS masses. 
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I. INTRODUCTION 



One of the main goals of the Large Hadron CoUider (LHC) is to search for new physics 
(NP) around or above the 1 TeV scale. Many new physics models predict heavy scalars 
carrying a color charge. Such scalars exist in supersymmetric theories [l|, |2|, Pati-Salam 
unification |3|, 1^, grand unified theories p-tZl, chiral color [sf, and top color [9J. Generi- 
cally, such particles can introduce unwanted flavor changing neutral currents (FCNC) and 
the usual expectation is that these particles must be quite heavy to avoid experimental 
constraints on FCNC. However, these constraints will depend on undetermined Yukawa 
couplings, and if suitable restrictions on these Yukawa couplings are imposed, the additional 
scalars can be surprisingly light. For example, FCNC constraints can be naturally avoided 
if one imposes Minimal Flavor Violation 10|, lUl on the Yukawa couplings of new physics to 
Standard Model (SM) fermions. Manohar-Wise 12| recently proposed an extended scalar 
sector with color-octet scalars (COS) that are also electroweak doublets, the unique repre- 
sentation consistent with MFV. ^ The existence of color-octet scalars of this type is weakly 
constrained by collider phenomenology because the COS couple most strongly to the third 
generation of quarks. Searches for new physics in final states with bbbb yield a rough con- 
straint of ms ^ 200 GeV, where ms is the COS mass, assuming that the COS Yukawa 



couplings to up- and down-type quarks are roughly equal [IJ]. Completely model indepen- 
dent constraints are even weaker. Ref. 15| concludes that masses of these particles could 
be as low as ~ 100 GeV and still be consistent with precision electroweak fits and collider 
data. 

In a recent paper 
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two of us argued that better constraints on the masses of COS 
can be obtained in searches for bound states of the COS. The COS can be pair-produced 
and have a strong attractive potential when they are in a color-singlet state. If the Yukawa 
couplings of the COS to SM fermions are 0(1) or smaller, this state can live long enough to 
form quarkonium-like bound states called octetonium. These bound states can then decay to 
pairs of electroweak bosons, e.g., 77, 7^'°, W~^W~ , etc. Thus the octetonium would appear 
as a resonance in these channels which have relatively small SM backgrounds. The couplings 
to gluons and electroweak bosons are fixed by gauge symmetry so the only free parameter 
in the calculation of the cross section is the COS mass. Ref. [l6[ calculated the production 
cross section for octetonia via gluon-gluon fusion as well as decay rates for a number of two- 
body decays to SM particles. A back of the envelope comparison of octetonium production 
via gluon-gluon fusion followed by decay to 77 suggested that the resonant cross section for 
this process would exceed the SM contribution for COS masses of 500 GeV or less. Thus, 
better constraints on the COS masses than those found in Refs. [M, ll5| could be obtained 
from null searches in these channels. 

The goal of this paper is to perform a more careful calculation of the process discussed 
in Ref. IGj by incorporating important QCD corrections that arise in the calculation of 
pairs of strongly interacting heavy particles near threshold. Many of the same issues arise in 



the calculation of tt, squark-anti-squark, and gluino pair production |l7H21|. There are two 



^ Other representations are possible if NP transforms non-trivially under the SM flavor group 13 1 



classes of corrections one needs to take into account. First, there are (partonic) threshold 
logarithms that appear in any production process characterized by a large partonic center- 
of-mass energy threshold including, e.g., Higgs production or Drell-Yan. In our previous 
the resummation of these logarithms for the production of a single COS was 
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performed using Soft-Collinear Effective Theory (SCET) |23l425l|. There we showed that the 
resummation increased the normalization of the total cross section by a factor of 2-3 for a 
COS with mass in the range 500 GeV-3 TeV. The additional effect that must be taken into 
account when two heavy colored particles are produced is the exchange of Coulomb gluons 
between the heavy particles in the final state. The exchanges scale as a^/f where v is the 
relative velocity of the heavy particles. In the threshold region, f ~ a^ graphs with Coulomb 
gluons must be resummed to all orders. The exchange of Coulomb gluons is responsible for 
the attractive potential between the COS when they are in a color-singlet state and gives 
rise to the resonant enhancement of the cross section when the invariant mass of the COS 
pair is close to that of the octetonium bound state. 

The outline of the paper is as follows. In Section II, we perform tree-level matching of 
the amplitude for gg — )■ S^S~ in the Manohar-Wise model onto SCET and Heavy Scalar 
Effective Theory (HSET). The resulting operator couples the COS, which are slowly moving 
and hence described by HSET fields, directly to the initial state gluons, which are described 
by SCET coUinear fields. In Section HI, we derive a factorization theorem for a{pp — )■ 
S^S^X). The cross section factors into a hard part (proportional to the square of the 
matching coefficient obtained in Section II), a soft function, and parton distribution functions 
(PDF's). Exchange of Coulomb gluons is included in the QCD Coulomb Green's function. 
In Section IV, we solve renormalization group equations (RGE) for each of the components 
in the factorization theorem. The resummed production cross section up to next-to leading 
logarithm (NLL) is obtained directly in momentum space using the methods of Ref . [26| . In 
section V, we extend our results to the cross section a{pp -^ S^S~ — )■ 77) and compare with 
the NLO SM calculation of a{pp — > 77) obtained using the program DIPHOX [27]. Before 
continuing, we wish to emphasize the universality of the factorization and resummation. All 
dependence on the model of NP is contained in the matching coefficients which enter the 
hard part of the cross section. The remaining steps of the calculation are independent of 
the model of NP. With suitable modification of the hard part, the results of this paper can 
be applied to any model of NP that contains COS. 



II. MATCHING S+S- PRODUCTION ON HSET/SCET OPERATORS 

At the LHC, the gg initial state gives the dominant contribution to pair production of 
COS. The leading order, O(a^), production processes are shown in the diagrams in Fig.[T]-(a). 
The couplings come from kinetic terms for the COS, 

Cs = -^S'^iD'rS' - ^mlS'^S^ (1) 

where D°f = 8^6°''^ + gf°'^'^A^ . We are interested in calculating the cross section in the 
threshold region, s ~ (2ms) ^, where s the momentum squared of incoming partons. In this 



region, the COS are moving slowly, but the initial state gluons are highly energetic. We 
integrate out the large scale ms and the match the diagrams in Fig. [I]-(a) onto operators 
with SCET collinear fields describing the initial state gluons and the HSET fields describing 
the slowly moving COS in the final state. After this matching, the production of S^ and 
S~ is accomplished via the operator depicted in Fig. [T]-(b) which is 

^™^ = £|^^'"'^''' + ^'"'^'"'^ {y^^n'Y {y^BQ\sryl)\s-yiy + o{a% (2) 

where n and n are lightcone vectors satisfying n^ = n^ = 0, n-n = 2, and -B^'^ is a leading n- 
collinear gluon field strength tensor, defined by B'^'^ = in^ g^^ G^^pJW^ = in^g'^yVl'^^G'^p^. 
^Til i^ related to B'^'^ by exchanging n and n. W^ is a collinear Wilson line in the adjoint 
representation 

Wf (x) = Pexpiig f dsn ■ A^(sn^)T'= j . (3) 

In Eq. (|2]), the heavy scalar fields are described by the HSET Lagrangian 

.chset = s:\v ■ iD.rs: - -^sT{Dirs:, (4) 

where v^ is the velocity and Ds is the covariant derivative including only the soft gluon field. 
The HSET Lagrangian can be obtained from Eq. ([1]) by making the substitution 

y/2ms V / 

dropping all terms in which the large phase does not cancel, and expanding to 0{l/ms). 
Note that we have decoupled soft gluons from the collinear and heavy fields by performing 



field redefinitions, so the soft Wilson lines 
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yf{x) = Fexplig f dsv ■ A^(sv'^)T^' j , v'^ = n^, n", v^". (6) 

appear in the operator in Eq. 02]) 

We find it useful to classify operators by the irreducible representation of color car- 
ried by the initial and final states. The possibilities are enumerated by applying 
8 (g) 8 = 1 © 8s © 8a © 10 © To © 27 to the initial and final states, and demand- 
ing that total color be conserved. There are eight operators that can contribute to 
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FIG. 1: Feynman diagrams for leading order color-octet pair production in full QCD (a) and 
effective theory (b). 



COS pair production via gluon-gluon fusion, and we will denote them by {Ri, Rf) = 
(1,1), (8s, 8s), (8s, 8a), (8a, 8s), (8a, 8a), (10, 10), (To, TO), and (27,27) [S], where Ri 
and Rf denote the irreducible representations of initial and final states, respectively. Using 
this operator basis the interaction Lagrangian in Eq. ([2]) is 



Ant = y^Cn,(/i)C>m(/i), 



(7) 



where each operator is given by 
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-if"'^'^, and we set A''^ = 3, where A''^ is a number of colors. 
All the color factors satisfy the orthonormality relation -E^birf-Eabcd ~ '^*"'- ^^ ^^®^ level the 
Wilson coefficients in Eq. (ITj) are {Ci,C2,Cs) = 7ras(6,6\/2, — 6\/3), and C3 = C4 = C5 = 
Ce = C7 = 0. In general, these color factor E^j^^^ ^' are defined to be J28| 
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(10) 



where C^^^ = {R,a\ab) are the Clebsch- Gordon coefficients between two different color 
spaces {R,a) and (a, 6). The -E^f,,,^^ vanish unless the two irreducible representations Ri 
and i?2 have the same dimension. 



III. FACTORIZATION FOR COS PAIR PRODUCTION 



For single COS production at threshold 221], the only degrees of freedom after integrating 
out the hard scale are the coUinear initial state partons, soft partons, and a single heavy 
COS. The interactions of the heavy COS with the initial state collinear partons via soft gluon 
exchange are equivalent to a time-like soft Wilson line. The resulting factorization theorem 
is a convolution of two parton distribution functions (PDFs) and soft function multiplied by 
hard Wilson coefficients. The factorization formula for COS pair production at threshold 



is similar, but Coulomb gluon exchanges between two COS in the final state must also be 
taken into account. These can be resummed to all orders using the QCD Coulomb Green's 
function, 



Gr(x,x',E) 



1 



Hr-E 
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C(x)C 
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(11) 



where Hji is the nonrelativistic Hamiltonian including the Coulomb potential in a specific 
irreducible representation, R, that for a pair of COS must be 1,8s, 8a, 10, 10, or 27. In 
the second equality of Eq. ( ITTj) . we have written the Coulomb Green's function in terms of 
the wavefunctions, ipl^{x.), of eigenstates with energy En. Our strategy for extracting the 
dependence of the cross section on the full Coulomb's Green function is to calculate the 
cross section for energy eigenstates and then use the second identity in Eq. (11 II) to infer the 
dependence on the full Couloumb's Green's function. We will also take into account the 
finite width of the COS by making the replacement E ^ E + iTs- 

The cross section in the threshold region for producing COS pairs is given by 
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where P^^ ^^^ ^^e incoming protons' momenta 



' Pn,n 



are the momenta of the partons, 



Mfc is a bound state mass, and the matrix elements Ail, T are defined to be Ai 



^m,k 



.{R,a) 



.{R,a) 



{0\. ' X\CmOm\PnPn) ■ Here the state \0\. ' ) is a COS pair in the color state (i?, a) and k 
refers to all other quantum numbers. The subscript m denotes the (i?j, Rf) quantum num- 
bers of the SCET operators and the summation is nonvanishing when the SCET operator's 
Rf is the same as the final state's R. For example, if we consider the final state with R = 8s, 
m can be either m = 2 (8s, 8s) or m = 5 (8a, 8s). 

The states X in Eq. ( 1T2|) consist of n(n)-collinear and soft partons, so the final state 
momentum and the phase space integral can be rewritten as px = PXn + PXrr + Pxs and 
Y2jf^ = X]x„ Sx- Sx ' respectively. The incoming parton momenta satisfy the relations, 
P{n,n) = P(n,n) — ViXn^x^)- Then the argument of the delta function in the last equahty of 
Eq. ( TT2l) becomes 
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where r] = 

2ms + E. 



Pxs^ Mk 



2ms + Ek, and the invariant mass of the COS pair is M = s 



1/2 



In order to derive the factorization formula in momentum space, we will insert into 0"/?^. 

1 = / dr]dyidy25{T] + ido)5 (yi - ^-p- j 6 iy2 - ^^— p- j , (14) 

where n ■ V{n ■ P) is a large label operator acting on n(n)-collinear fields, and the partial 
derivative, ido, gives the energy of soft partons. Using the definition of Oi in Eq. ([H]) and 
the completeness relation |X)(X| = 1, we write aptf as 
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where s = yiy2S. In the third line of Eq. ( !T5|) . we introduced the following notation 

n-V 
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(17) 



The parton distribution function (PDF) for the gluon in Eq. flT6|) is defined by 

{Pn\B^'''B^^'[y]\Pn) = aTS'^'^^^^fMy)- 

The same equation, with n and n exchanged, defines for the n-collinear gluon PDF 
In Eq. ( !T6|) . the bound states \0n ) are defined in terms of the COS states by 



(18) 
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C(k)|^r(k)5-^(-k)), 



(19) 



where ip^O^) ^^^ the nonrelativistic wave functions in momentum space. The COS single- 
particle states annihilated by the HSET field, S*", are related to the states annihilated by the 
field S"" by \S!^) = {l/\/2ms)\S°'). In case of the final state with two identical particles such 
as \S^S^), the right side of Eq. flT^ should be divided by a/2. Finally the matrix elements 
for the color-octet scalars in Eq. ( !T6|l can be written in terms of the wavefunctions at the 
origin using 
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Here we identified M ~ Mk in the second equality, ignoring 0{1/M) corrections. Using the 
relation 



inserting Eqs. ( ITTl) . ( TT31) . and ( 1191) into Eq. ( TlGl) . and using 

5^%^ - M,2)|^f^(0)p = ^lmG^/0,0,i?) , (22) 

k 

we see that the cross section is proportional to the imaginary part of the Coulomb Green's 
function evaluated at x = x' = 0. Because the COS are unstable, we make the substitution 
E ^ E + iVs- 

Combining Eqs. (Il6]), ([20]), and (ES]), we then obtain 

(jR^iw^SSX) = ^^e ^j^2 _ iy J2 J dr]dyidy2s\CR^.R^{M, ^j)\'' 

X \mGR^.{(),Q,E + iTs)fg/p{yi)f,/p{y2)SR^,R^{7]), (23) 

where we have slightly modified our notation by replacing Ylm "with YIr ^^"^ ^m with 
Cr.^Rj,. The function Sr^^r^, is defined by 

Sr^Mv) = ^dM E^f^>Ei;;^^^Ei^f^> (24) 



X (0 
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When we introduce the variable z = M'^/s, which goes to 1 at threshold, the soft mo- 
mentum 1] can be rewritten as 

£l/2 

ri = s^/^-M = s^/\l-z^/^)r^^il-z), z^l. (25) 

Using the relation yiy2 = t/z and replacing j df] ^^ —J dM, we find that the differential 
scattering cross section is 

^(pp^S+S-X) = J2HR,,Rf{M,l^F)j^-r^ lmGR^.iO,0,E + tTs,lXF) (26) 
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where ixp is the factorization scale. Fir j z) is a convolution of two PDF's, 

F{x,^f)= \ —fg/piy,i^F)fg/pix/y,fXF), (27) 

and the dimensionless soft functions, Sr-^Rj:{1 — z), are defined to be Sr.^Rj_.{1 — z) = 
{s^l"^ /2)SR^^R^{ri), so that Sr-^r^{1 — z) = 5{1 — z) at tree level. Finally, the hard func- 
tion Hr^^Rj is 

gi?.,i?/M,/i) = 16 /^2_i)2 • (28) 



This factorization formula is one of our main results, and can be extended to other processes 
with different initial states such as qq and qq. Note that to obtain the scattering cross section 
for the production of two identical particles, such as S^S^ , the cross section should be divided 
by 2. If we restrict the sum over bound states to the ground state of the singlet channel 
[Rf = 1 and k = corresponding to the state with principal quantum number n = 1 and 
/ = 0), and use the tree-level soft function, Sr-^ji^ = 6{1 — z), Eq. (12^ becomes 

4%P -* 01) = ^^f^'ll°l, Mm'rF(ry (29) 



This reproduces the tree-level cross section for pp — )■ O^ in Ref. [16 . 

In the scattering cross section, the hard function, H, the soft function, S, and the 
Coulomb Green's function in Eq. ( 126|) should be evaluated at renormalization scales la- 
belled fiH, /^s, and He, respectively. These scales are chosen so that large logarithms are 
minimized. Large logarithms are resummed by evolving the hard function from fin to iip, 
the soft function from fis to ^p, and the Green's function fj,c to fj,F- This is described in 
the next section. 



IV. RESUMMATION AND NUMERICAL RESULTS 

In this section we calculate the resummed scattering cross section to NLL accuracy with 
leading order (LO) Wilson coefficients. This approximation, called NLL+LO, includes all 
0{1) terms when the large logarithms are counted as an inverse power of a^, so corrections 
to NLL+LO are 0(0;^) suppressed. Next-to-leadin g or der (NLO) Wilson coefficients have 



not been calculated for COS pair production. Refs. [29|, ISOj have observed that 7r^-enhanced 
NLO contributions, which can be inferred from the imaginary parts of anomalous dimensions 
as we will see below, are numerically similar in size to the complete NLO «« correction. So 
we will include this contribution in our numerical results, which, based on expectations from 



previous calculations of Higgs production [30|, should provide a result numerically consistent 
with a full NLL-I-NLO calculation. 

As shown in Sec. [TTl the only nonzero LO Wilson coefficients for gg — )■ 5* 5" are Ci, 
C2, and Cs, which correspond to the initial and final states (1,1), (8s, 8s), and (27,27), 
respectively. Computing anomalous dimensions for the NLL resummation in each of these 
channels is straightforward, and the results are 



a. 



iM,) - -(s^o^- (s)'rO'"3]i^-S< (») 



where Ca = N^, B^ = 2/3o, S^ = 2Ca + 2/3o, B^^ = 16 + 2/3o, and /3o is the first coefficient 
of the QCD beta function. Here, F^ and F^ are the first and second coefficients of the 



cusp anomalous dimension of Wilson lines in the adjoint representation: Tq = ACa and 
r^ = 8A''c[(67/18 — ■n'^ /Q)Nc — 5?2//9], where nj is a number of flavors. 

From the anomalous dimensions we can infer the form of the double logarithms in the 
Wilson coefficients, which are 
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These lead to large vr^-enhanced corrections when evaluated at the scale fi = M 



~ic{;;,3^(M)pexp(^c4i,i,-i; 



3. 
In the second line we have exponentiated, the Tr^-enhanced terms. This is a consequence of 



evolving the renormalization scale to a complex value so as to minimize the logarithms [22 . 
Interestingly, the Wilson coefficient for the 27 channel is suppressed when the vr^-enhanced 
contribution is included. 

In the resummed cross section, we use the tree level values for the soft functions in 
Eqs. (!24l) and ( !26l) . However, we need to evolve the soft functions from the soft scale, fis, 
to the factorization scale, fip, and to determine the appropriate fis we will use the one- loop 
expressions for the soft-functions: 

Si,i(l -z,fi) = 5(1 -z) + —N,A{1 -z,ii), (35) 

vr 

S8s,8s(l-^,/i) = 5(l-^) + giV,(2A(l-^,/i)+5(l-z,/i)), (36) 

^27,27(1 -^,/i) = 5{l-z) + ^N,{2A{l~z,ii) + ^-B{l-z,ii)), (37) 

where the coefficient function A{1 — z,fi) is obtained from soft gluon exchanges between 
3^„ and 34: or yn{n) and y^, and 5(1 — z, /i) from soft interactions between 3^^'s. These 
coefficient functions are 

'v^^-^Va-.)-21nii— i- + 4r^ 



A(l-z,a) = ll\n' Jl- - —]5(1 - z) -2\n^- — ^^ + 4 ^^^^^ ^ , (38) 

^ '^^ V2 M2 4/ ^ ' M^l-z)+ \ 1-z J+ ^ ' 

5(1 -^,/i) = (1^^ + 2)^(1 -^)-(r^' (39) 

where the standard plus distributions are used, and ultraviolet (UV) poles have been ab- 
sorbed into counterterms. Note that these expressions are infrared (IR) finite. The general 
form for the NLO soft function for the process /1/2 — )■ F, where Ji and I2 denote the color 
representations of the initial partons and F denotes the irreducible representation of the 
final two heavy particle states, is given by 

:5(i„i,),f(1 -z,fi) = 6{1 - ^) + ^ ((Ci, + CiJA(l -z,fi) + CfB{1 - z, /x)) , (40) 
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where Ci^ ^ and Cp are the quadratic Casimir operators for the initial and final represen- 
tations. Our result agrees with Ref. 
coordinate space. 



28| . where the computation has been performed in 



The Coulomb Green's functions in Eqs. fl22|) and fl26|) are 17|, |31 



GRAO,0,E + iTs,fi) 



Arc 



CRj^ml 



Ik \2msvJ 2 ^^ ^ 



(41) 



where k, v, and ip are 



2^ ' 



E + tTi 



ms 



^{z) 



1E + — lnr(2;). 
dz 



(42) 



Here E = M — 2ms, 1e is the Euler gamma, and T{z) is the Gamma function. The 



-a.C 



Rj/r, where Rf 



Crj are the coefficients in the LO Coulomb potential, Vc,Rf{r) - 
refers the representation of the COS pair. For COS pairs, Ci = Nc, C^s = Nc/2, and 
C27 = —1, so the COS pairs in the 1 and 8s feel an attractive force while COS pairs in 
the 27 feel a repulsive force. The appropriate scale for the Coulomb's Green's function is 
/ic* ~ ^^5^ ~ '"^5C*%tts(/^c)) where v is the relative velocity of the COS. In the resummed 
cross section the Coulomb Green's function needs to be evolved from the scale nc to the scale 
Hf, as indicated in Eq. (l26i) . However, the Coulomb Green's function anomalous dimension 
starts at O^af) so its evolution can be neglected in a NLL calculation. 

The renormalization group equations (RGEs) for the hard functions, soft functions, and 
PDF's are solved directly in momentum space using the methods of Ref. 26|, |32|, |33|. The 
details of the calculation are very similar to the calculation of the resummed cross section 
for single COS production in Ref. [22| so we simply quote our result for the differential cross 
section: 



dcTR^. 
~dM 



{pp -^ s+s-x) 



M 



{2ms 



lmGR^.{0,0,E + tTs,fic) 



f dz 

XT —Vr^{z,M,hf)E{t/z,hf), 

J z Z 



(43) 



where the resummation function, Vr^.{z,M, ^xf), is given by 



VR^.(z,M,fif) 



Ri 



-2'yEV 



U-2r, 



n^v) 



(44) 

Here SR.^Rj,{dr,, fis) are the Laplace transforms of the soft functions, and the evolution func- 
tions UR.^R^.{fiH, fJ'S, fJ'p) are multiplicative factors that come from evolving the hard func- 
tions from the scale fiF to the scale fin and the soft functions from the scale fiF to the 
scale fis- Up to NLL accuracy the representations of the initial and final states are the 
same, so below we will simplify our notation by replacing /r-^r^ with /r^ where / represents 
either a hard function, soft function, or evolution function, and suppress the summation 
over Ri in Eq. (jSj). For the NLL resummation, the auxiliary parameter rj is defined to be 
1] = {T^//3o)ln{as{fif)/as{fJ.s)) as in Ref. [26 . 
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FIG. 2: Mass distribution of the scattering cross section for pp — ?■ S~^S~X near threshold for ^/s = 
7 TeV. Upper (Lower) fihed regions with blue (Yellow) color represent da-i/dM {da^^/dM) with 
the soft scale varied between /x^/^ < l^s ^ l^^s- 



The NLL expressions for URj,{fiH, fJ^S: fJ^p) are 



B 



InlJji {ij.h,IJ'S,I^f) = In 4S'?7nll (/"//, /U5) + -5— In 



Rf , Osif^s) . Bg asifJ^F 



+ ^ln' 



(45) 



/3o asifin) /3o asifJ-s) 
where Bg = 2(3q and B^ are defined in Eqs. fl30|) . fl3T|) . and (|32|) . The function 5'[/nll(/^15 /^2) 



IS 



^f/NLL(/.i,/.2) = ^[^(l---lnr) + (^-^)(l-r + lnr) + Ain^ 
4/3^ La^(/ii) V r / Vr^^ /3o/ 2/3o 



(46) 



where r = as(/i2)/as(/^i)- 

We will choose the hard scale to be fin = M, and use the second line of Eq. ( l34l) for 
the Wilson coefficients in the hard functions so that the large vr^-enhanced contribution is 
included. To resum logarithms of 1 — 2;, the soft scale should be set to /X5 = M(l — z). 
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FIG. 3: Mass distribution of the scattering cross section for pp — t- S~^S X near threshold for \/i = 
14 TeV. Upper (Lower) fihed regions with blue (yellow) color represent da\/dM (dagg/dM) with 



the soft scale varied between fj,g < /X5 < ^^ 



However this choice gives divergences in the z integral since the running coupling will cross 
the Landau pole as z — )■ 1. Instead we chose the scale ^s so that the higher order corrections 
to the soft function are perturbatively small. In order to do this, we define two soft scales, 
lig and fig^. The scale fi^^ is defined by starting from /i^ = fin and lowering fis until the 
0(as) correction is less than 15%. The scale /i^^ is chosen so that the one-loop correction is 



30,13^ 



minimized. The soft scale fis is then defined to be the mean of fig and fXg^ 

In the Manohar-Wise model [l2|, the dominant decay modes for S*^ with masses greater 
than 200 GeV are S+ -> tb {S' -^ ib), so the width of the COS is gven by 



r^ ^ r(^+ -^ tb) 



\vum 



tb\ 



167rm| 



V V 



2\2 



{rris-m^) 



(47) 



where v is the vacuum expectation value of Higgs and rju is an unconstrained complex 
Yukawa coupling. Eq. (H71) is used to determine F^ in the evaluation of the Coulomb Green's 
function. When \riu\ is smaller than 1, the COS live long enough to form bound states called 
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octetonium jl6|. Figs. [2] and [3] show the differential cross section as a function of the invariant 



mass of the COS pair. For the numerical results we employed NLO CTEQ5 PDF set [34 
Fig. [2] shows the cross section at y/s = 7 TeV for m^ = 350 GeV and 500 GeV, and for 
two values of the Yukawa coupling, \riij\ = 0.5 and 1.0. Fig. [3] shows the same for a/s = 14 
TeV. As seen in Figs. |2] and El the octetonium appears as a resonance 10-15 GeV below 2ms 
that is clearly visible in the 1 channel. In the 8s channel, there is a small peak in the cross 
section just a few GeV below 2ms. This peak is so broad for \riu\ = 1.0 that it is barely 
noticeable, but the peak is visible when the Yuklawa coupling is \riu\ = 0.5. The scattering 
cross sections in the 27 channel do not have peaks and are negligible compared to the 1 and 
8s channels, so we have not included them in Figs. [2] and Fig. [31 



V. SCATTERING CROSS SECTION FOR pp -^ S+S^ -^ 77 



Ref. [IGJ argued that the process pp — )■ SS — )■ AB, where AB represents a pair of SM 
electroweak gauge bosons, e.g., W^W~ , Z'^ Z^ ,'-/'-/ , or 7^°, are promising channels in which 
to search for octetonium. Near the vicinity of the octetonium resonance there is a peak 
in the cross section which can exceed the SM background for these final states. This is 
in contrast with final states like gg or ti where we expect the QCD background to greatly 



exceed any signal from octetonium. In Ref. 16|], a simple estimate for the cross section for 



pp — )■ S^S*^ — )■ 77 in the vicinity of the octetonium resonance was compared with the SM 
background. For octetonium with mass ^ 1 TeV (m^ < 500 GeV) the cross section at the 
LHC at \/i = 14 TeV was found to exceed the SM background for this process. Therefore, 
searches for 77 resonances could either reveal these novel heavy states or provide much 
better constraints on the allowed masses of COS, which are currently only constrained to 
be > 100 GeV |15i]. The point of this section of the paper is to improve upon the results 



of Ref. [16[ by performing a resummed calculation of the invariant mass spectrum for the 
photons produced in pp -^ S^S^ — )■ 77 in the vicinity of the octetonium resonance, which 
is compared with the SM prediction for the 77 invariant mass distribution. 

Below the threshold 2m,s (and ignoring the widths of the bound states), the cross section 
for pp — )■ S~^S~ — )■ 77 can be written as a sum over contributions from individual Ol states, 

a,{pp ^SS^ 77) = E^i.'^(PP ^ ^'^^^ X'^^^n^jtv (^8) 

Note that only color-singlet resonances can decay to the final state 77. Using the factor- 
ization formulae in Eq. ( H3l) . integrating over M, and writing ImGij^=i(0, 0, E) as a sum of 
^-functions times wavefunctions squared, as in Eq. ( p2|) . we can write cr^^kipV ~^ 0\X) as 

<^i,k{pp ^ OlX) = ^^-— ^ j dMH^{M, ^i)TM' (49) 

x6iM'-M'MliO,^^)\'J^'^S,il-z,^^)F[-^,^?|. 
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The LO decay rates for Ol -^ 77 are [16 



rl{ol-^^j) = ^^\ijlm'- (50) 

Finally, we must allow for a finite width for each of the bound states, O^. We do this 
by replacing 6{M^ - M|) with the Breit-Wigner {MkTl^^^jTT)/{{M^ - Mlf + M|(r^ 



\2\ 



toty 



Then we can simplify Eq. fH8|) with the substitution 



r£ot \^m? \^im 



t^o ^ 4M7r E-Ek + tTl^j2 E - Ek - tTltoJ^ 



4Mn^^E-Ek + iTl,j2 E - E, - zTl,j2 
^ %^\Gi{0AE + zTlJ2)\'. (51) 

In the second line we replaced Tl ^^^ with Fq ^^^ so we could write the final result in terms 
of the Green's function. The corrections to this approximation are not important near the 
resonance of interest, and small except near the other poles of the Green's function which 
should not be important for our calculation. Then we combine Eqs. (1491) . (1501) . and (ISTIl to 
obtain 

ai(pp^5+5-^77) = ^^y"^i/i(M,/i)|Gi(0,0,E + zFg„,/2)p 

xry"'^:Si(l-z,/i)F(^,/i). (52) 

Using Eq. fl52|) . we compare the cross section for pp — >■ S~^S~ — )■ 77 to the SM background 
pp — )■ 77. The cross section is computed in the vicinity of the 0^{= Oq) resonance in the 
Manohar-Wise model. In this model, the width of O^ depends on a scalar coupling, Ai, which 
appears in the coupling of a COS pair to the SM Higgs boson jl2|. We have set this parameter 
to Ai = 1. Explicit expressions for the decay rates for O^ — t- gg, ti, W~^W~ , Z^Z^, 77, and 



hh can be found in Ref. UM, and these have been used to calculate the total width, Fj^^^. 



In Fig. m we compare the 77 invariant mass distribution near the peak of the resonance O^ 
with SM backgrounds for pp — )• 77. The SM background cross section, dasu/dM, is the sum 
of NLO calculations of qq — )■ 77 and gg — )■ 77 with a rapidity cuts of |?7i,2| < 2.4. The K- 



factor has been computed using the program DIPHOX [27[. The cross section is computed 
for v^ = 7 TeV with COS masses of 250, 300, and 350 GeV, and for v^ = 14 TeV with 
COS masses of 300, 400, and 500 GeV. We see that the resonant cross section exceeds the 
SM contribution when ms < 500 (350) GeV for ^/s = 14 (7) TeV, confirming the conclusions 
of Ref. [16]. We have used 'qu = 1 in our calculation, for smaller t]i/ the resonance peak 
is more narrow and visible. Note that the octetonium, O^, which is composed of a pair of 
electrically neutral COS, is significantly narrower than O^ when \r]ij\ = 1.0 and rris < 700 
GeV. Therefore, this should appear as a narrower resonance in channels into which it can 
decay, such as W~^W~ and Z^Z*^. It would interesting to extend the results of this paper to 
other final states with electroweak bosons. 
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FIG. 4: Mass distribution of the scattering cross section da\/dM{pp — ;■ S^S~ — ;■ 77) near the 
resonance O^ (= Og) versus the Standard Model background dasM/dM{pp — )• 77) for (a) ^/s = 
7 TeV and (b) y/s = U TeV. 



VI. CONCLUSIONS 

The LHC will explore physics beyond the TeV scale. One possibility for new physics that 
may be discovered at the LHC is the existence of heavy color-octet scalars (COS). In this 
work we have extended our previous analysis of single COS production [2^ and considered 
the production cross section of two COS which bind together through Coulomb interactions 
to form a bound state called octetonium [16[ . This bound state can decay into two photons, 
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providing a resonant signal above the SM di-photon production cross section. We established 
a factorization theorem for this production process using SCET and HSET, then performed 
a next-to-leading logarithmic partonic threshold resummation directly in momentum space. 
Our factorized cross section is independent of the specifics of the underlying NP theory 
responsible for the production of COS. In this paper, we focused on the Manohar-Wise 
model of COS, but the calculation can be easily extended to pair production of heavy colored 



particles in other models, e.g. stoponium in supersymmetry [35H37J or pairs of Kaluza-Klein 
excitations of quarks and gluons in models of extra dimensions, by a suitable modification 
of the matching coefficient at the high scale, 2Ms- We find that the resonant cross section 
exceeds the SM contribution when m^ < 500 (350) GeV for -^ = 14 (7) TeV. Searches for 
di-photon resonances at the LHC will either discover COS resonances or greatly improve 
existing bounds on COS masses. 
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